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We study the survival probability and the first-passage time distribution for a Brownian motion 
in a planar wedge with infinite absorbing edges. We generalize existing results obtained for wedge 
angles of the form n/n with n a positive integer to arbitrary angles, which in particular cover the 
case of obtuse angles. We give explicit and simple expressions of the survival probability and the 
first-passage time distribution in which the difference between an arbitrary angle and a submultiple 
of 7r is contained in three additional terms. As an application, we obtain the short time development 
of the survival probability in a wedge of arbitrary angle. 

PACS numbers: 


I. INTRODUCTION 

The survival probability, which is the probability not 
to have reached a target up to a given time, is a key 
observable in the study of Brownian motion. It is the 
cumulative of the first-passage time distribution, i.e. the 
distribution of the time to reach a target. This quantity 
[TH5! , and as a first step the mean first-passage time [B] , 
is a standard way to quantify the efficiency of a search 
process. It is for example also involved in the calculation 
of the covered territory 13 la and the mean perimeter of 
the convex hull of a Brownian motion [!j] . 

Determining the survival probability and the first- 
passage time distribution is of particular importance in 
the wedge geometry. Recent related studies include last- 
passage times determination m , extensions to anoma¬ 
lous diffusion |TT] (and in particular to fractional Brow¬ 
nian motion m) and applications to virus trafficking in 
cells m- Interest in this geometry resides in part in the 
possibility to map one-dimensional diffusion-controlled 
reaction processes on a wedge domain POH5]. One 
important example of this mapping is the Fisher-Gelfand 
three-particle problem poms]- Consider three diffusing 
particles on a line, with diffusion constants D\, D 2 and 
D 3 . Given the starting positions x±, aq and x 3 , what 
is the survival probability of the middle particle up to 
time t, that is to say the probability that it has not 
met the two other particles up to time f? By writing 
the Fokker-Planck equation in coordinates yi = x\ — X 2 
and IJ 2 = X 2 — X 3 , and after some transformations, this 
problem reduces precisely to the problem of the survival 
probability of a single Brownian motion in a 2D wedge 
with top angle a = 2 arctan[y /(1 — q)/(l + 7 )] where 

7 = D 2 /\J (Di + D 2 )(.D2 + D 3 ). 

In the general case, the survival probability for regular 
diffusion in a wedge domain is written as an infinite sum 
of special functions. In the past, special attention has 
been devoted to the analysis of the large time behavior 
of the survival probability, which displays a power law 
decay with an exponent continuously depending on the 


wedge angle p. On the other hand, the analysis of the 
short time behavior seems to have been left aside. 

Recently, and in contrast with the standard expression 
involving an infinite sum, compact analytical expressions 
of the survival probability and the first-passage time dis¬ 
tribution have been obtained for special values of the 
wedge angle by using the method of images [TB]. These 
expressions have been extended to biased diffusion 113. 
However, these results have been limited to specific wedge 
angles of the form ir/n where n is a positive integer. In 
particular, they do not apply to obtuse angles. 

Here, focusing on unbiased diffusion, (i) we general¬ 
ize these results to arbitrary wedge angles; Note that 
this covers in particular the case of obtuse angles, which 
has proven to be an essential ingredient in the context 
of diffusion growth processes P jTSJ; Beyond this, it is 
crucial to know the survival probability in a wedge of 
arbitrary angle to solve the Fisher-Gelfand problem pre¬ 
sented above for general diffusion coefficients; (ii) The 
expressions presented here take the form of a finite sum 
over generalized images plus an integral only involving 
elementary functions, which vanishes for wedge angles of 
the form 7 t /(2 p +1); This structure thus explicitly un¬ 
derlines the difference between a wedge angle rr/n and 
an arbitrary one; (iii) As an application, we show that 
the short time behavior of the survival probability is con¬ 
veniently obtained from this expression. 

The manuscript is organized as follows. In Sec. II, we 
remind the standard expression of the survival probabil¬ 
ity and give the main results of this paper. In Sec. Ill, 
we provide the derivation of the alternative expressions 
of the survival probability and the first-passage time dis¬ 
tribution for any wedge angle - the most technical part of 
the derivation appears in appendix - and compare them 
with the existing results for special n/n wedge angles. In 
Sec. IV, we give the short time asymptotic development 
of the survival probability. Finally, in Sec. V, we draw 
our conclusions. 
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II. BASIC EQUATIONS AND MAIN RESULTS 

Let us consider a Brownian motion in a planar wedge 
of top angle a with two infinite absorbing edges, start¬ 
ing from a point (ro,</?o) (see Fig. [l]). The propagator 



Figure 1: Geometry of the wedge. The starting point is 
(ro,<po) at time 0 and the position of the walker at time t 
is (r,p). 


P(r, p, t\ro, po) satisfies the diffusion equation 


dP 

— = DAP = D 
dt 


d 2 P 1 dP 1 d 2 P 
dr 2 ^ r dr r 2 dtp 2 


(1) 


with the following initial condition and boundary condi¬ 
tions: 

P(r, p, 0|r o , p 0 ) = S(r - r 0 ) = —S(r - r 0 )5(p - p 0 ) 

r o 

P(r, 0, t\r 0 , ip 0 ) = P(r, a , t\r 0 , p 0 ) = 0 ( 2 ) 

The exact solution of this problem is known mm 

+oo 

. /mrip\ . (mript 

5in I - 1 sin I - 

a ) 


1 -l-oo 

P(r, tp,t\ r o, po) = ^E sin (^) sin ( 


Xli 


/ r 0 r \ 
V2 Dt) 


exp - 


4 Dt 


( 3 ) 


The survival probability, defined by 

/*+oo pOi 

S(t\r 0 ,<p 0 )= dr r dp P(t,r,p\r 0 ,Po), (4) 

Jo Jo 

can be computed by integration by parts, as a function 
of a rescaled variable y = r\j ( 8 Dt) and the initial angle 
Po 


s(y,p o) 


2 M e -y Y' sin (( 2w + 1 )pr) 

V 7T 2m +1 

—n 


x \\. v (y) +I^+i (y)] 


( 5 ) 


with v = ( 2 m + l)n/( 2 a) — 1 / 2 . 

While this expression is well-suited to analyze the 

large time (small y) behavior of the survival probability 

S(y,ipo) oc y n A 2a \ it is difficult to extract the small 
y-t o 

time (large y) behavior. Indeed, using bluntly the be¬ 
havior of the modified Bessel function to large argument 

i Ay) ~ -f=(i-^^) ( 6 ) 

y->oo y/ 2 ny \ 8 y ) 


leads to S(y,<p o) —> 1 but does not allow to obtain 

y—> oo 

higher order corrections. The main purpose of this pa¬ 
per is to provide this small time (large y) asymptotics. 
We proceed to provide a compact alternative expression 
of the survival probability (summarized in Eqs. (25) and 
(26)) that is more suitable for the large y analysis. The 


result for the large y asymptotics are summarized in 
Eqs. (391, (401 and (41). In particular, we show that 


for large y, 


s(y,<p o) ~ i- 

2 /—>■ oo 


erfc (i/^/nun (p 0 , (V 


and also provide further subleading corrections. 


III. DERIVATION OF THE SURVIVAL 
PROBABILITY AND THE FIRST-PASSAGE 
TIME DISTRIBUTION 

The starting point to establish this alternative expres¬ 
sion of the survival probability is the integral representa¬ 
tion of the modified Bessel function I„ 

I v(y) = — [ dip e ycos>p cos(vp) 


sin(i/7r) 


r+oo 

/ d M e- ycosh “-"“. ( 8 ) 

Jo 


Plugging this form into Eq. (Js]) , the survival probability 
becomes 

/9 \ 3 / 2 

S(v,<Po)=[^) Vye~ y (Ai+A 2 ) (9) 


with 


/*7T +00 

A\= dp e ycosv (cos(i^y) + cos((^ + l)y>)) 
"'° m= o 


sin (( 2 m + 1 )^/ 


2m + 1 


P - 


= 2 dpe vcosv cos y -B 1 
Jo 2 


( 10 ) 


with 


B 1 

+oo 

= 5Z cos ( 

m—0 


irp\ sin (( 2 m+ 1 )^) 


2 a J 


2m + 1 


( 11 ) 


and, replacing v with its expression, 

p+oo 


r-r oo 

A 2 = due~ vcoshu (e~ u -l) 
Jo 


/ 0 

+oo 


V— . , sin (( 2 m + 1 )^) _ vu 
y j sm(i/ 7 r)-^- a e 


m —0 


2m + 1 


P+OO 


= 2 


du e _y cosh “ sinh —B 2 ( 12 ) 
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with 


. al 3a 

(2k+l)a 

+00 , 2 \ 

B 2 = £cos (2m + l}|- 

m—0 ' J 

sin((2m+l)^) e _ (2rn+1)s 

2a 

2ka (2k+2)a 

2m + 1 

(13) 

a| 3a 

(2k+l)a 



2a 

2ka 1 n (2k+2)a 


A. Calculation of the term Ai 

We show in this section that the term Ay can be ex¬ 
plicitly calculated. We first note that the sum By can be 
rewritten as 


Figure 2: The integer k is defined in Eq. ( |16| ). The number 
(2k + 2)a can be either smaller than n (first case), or greater 
(second case). In the first case, the interval [(2k + l)a, 7r] has 
to be cut again into the two intervals [(2k + l)a, (2k + 2)a] 
and [(2fc + 2)a, n]. 


+oo 


Bi= da/ cos((2 m + l)a/) cos ((2m + 1)^) 


m—0 


(14) 

with if varying from 0 to n. Then, we make use of the 
following formula 


+oo J- 

Y cos ((2m + 1)cos ((2m + = ^(y - z), 

m—0 

(15) 

valid for 0 ^ y ^ L/2 and 0 ^ z ^ L/2, for L = n. Note 
that special attention must be paid to the ranges of x' 
and nip/(2a) in Eq. (141 in order to use Eq. ( p~5| ). The 
condition 0 ^ x' ^ 7r/2is always respected as ifo ^ a/2. 
In turn, the respect of the condition 0 ^ irif/(2a) ^ tt/2 
for if in [0,7r] depends on the value of a. 

We first address in detail the case a ^ n. The key 
point is to cut the interval of variation of if into well- 
chosen intervals on which, after periodicity and parity 
manipulations on the cosine, the formula (151 applies. 
We define the integer k such that 


can then write the general relation 

[ dip f(ip) = f dtp /(VO 

Jo Jo 


Cl 


k 

+ V / d ^’ (/( 2 J“ - VO + /( 2 ja + VO) 

3=1 J ° 


+ 


' max(0,(2fc+2)o;—7r) 


c 2 =c^+c+ 

d'lp f((2k + 2)a — VO 


/•max(0,7r—(2fe+2)ci!) 

+ / dif) f((2k + 2)a + V') ■ 

Jo 


(17) 


with C g" that equals zero in the case where 
(2k + 2)a ^ 7 r. 


(2k + l)a ^7r < (2k + 3)a, (16) 

that is to say k = \jr/(2a) — 1/2J. We cut the interval 
[0,7r] into three: 

• the part [0, a], 

• the part [a, (2k + l)a], 


• the part [(2k + l)a,7r]. 


Referring to Fig. [2] we can see that two cases arise, de¬ 
pending on whether (2k + 2 )a is smaller or greater than 
7 T. If (2k + 2 )a < 7r, we cut again the last interval into 
two parts [(2k + 1 )a, (2k + 2)a] and [(2k + 2)a, 7r]. We 


The calculation of the integrals C\, C 2 , C 3 and C 3 , 
carried out in appendix, leads to the final expression for 

A! 

. / 7 r \ 3 / 2 e y 

k 

erf ^v/2ysin^ 0 ) + 5^(-l) J erf (y^ sin(ja + y> 0 )) 
7 — 1 

erf [y/2y sin(ja - ip 0 fj 


+(-l) fe+1 erf (y/2y sin (min ((ft + 1 )a + ip 0 , |))) 


— erf ( \J2y sin (min ((k + 1 )a — ipo, 


(18) 


In the case a ^ 7r, the term Ay is obtained by following 
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the same lines 

Ai = (|) *A__ erf (\/2y sin (min (y 0 , |))) • (19) 

B. Calculation of the term A 2 

The sum B 2 can be rewritten as 


Bo = - 


where 



-BE 

/ dx' (V) + 

/ dx' Co(x') 

Jo 

/ 0 


( 20 ) 


+OO 


ciix') = ^2 cos ( (2m + 1) f a;' ± — j j e < 2m+1 )5a 

ra— 0 ' ' ' ' 

( 21 ) 

After simple algebra, we get 


Cf(x') = Re 


e (2m+i)(*(*'±§|)-f“) 

(*'±fi) 


_m— 0 

2 COS 


2 sinh (^~) ^ sin 2 fx'± ^ 


( 22 ) 


We integrate over x' 


Bo =- 


arctan 


■ arctan 


sinh 2 (g) 


, sinh (e) ) 

' sin (S(^ o-f)) 


sinh(^) 


(23) 


and finally, 


1 r + °° 11 

Ao = - du e _1/cosh “ sinh — x 

“ 2 Jo 2 


arctan( ™(f ^ + f» ) + arctanf (I ~ Ilf 

v “ h <e) ) \ siab (s) 


in an acute wedge of top angle a ^ 7r 
S(y,Vo) = erf sin 

k 

+ ^(-l) J erf (-\/2y sin (ja + <po)^ 
3=1 


- erf yy2y sin (ja - ip 0 ) 

+(—l) fe+1 erf (-y/2ysin (min ((fc + l)a + ipo, —))) 

- erf (\/2y sin (min ((fc + l)a - <p 0 , 0)) 

- 5 * 

arctan ( ”"(% + ?)) ) + arctan Mi ' «>' 




sinh(l) 


(25) 


and in an obtuse wedge of top angle a ^ tt 
S( y,<Po) = erf (v^sin (min (y 0 , |))) 

+ (^) ' y/y-^- J+°° du e _!/cosh “ sinh ^ x 
arctanf sin (i (<*> + *)) ')+ arcta nf 

V sinh (e) ) V sinh (e) 


(26) 


with k = L7r/(2a) —1/2J and y = r^/fiDt). We condense 
Eqs. (|25| and (261 into the following expression 


S(y,<Po) = erf (\/2y sin (min (y 0 , ^))) 
fc +1 

+ ^(-i) J erf (v/2ysin (min (ja +y 0 , ^))) 
erf (v/2y sin (min (ja - y 0 , y)) 


3 = 1 


+ 


3/2 




=-v r+°° 


du e _y cosh “ sinh — x 
2 


arctanf 3in <° ( ? + j» ) + arctanf Si " ( ° (,| j-f» 

V “‘‘(S) ) V sinll (S) 


(27) 


(24) 


C. Survival probability 

1. Expression for an arbitrary angle 

We gather the results of the two previous subsections 
and give the final expression of the survival probability 


This expression is valid for any arbitrary wedge angle 
and generalizes the expression of Dy and Esguerra that 
requires wedge angles of the form 7r jn with n an integer 
m, as we proceed to show in the next paragraph. Note 
that the case of obtuse wedges was not covered by their 
approach. 


2. Particular cases n/n 

We check that in the particular cases of wedge angles 
of the form n/n with n an integer, we recover the expres- 
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sions of Dy and Esguerra mm- First, if n = 2p + 1, 
the integral part of the survival probability disappears, 
as easily seen in Eq. In this case, k = p, and 

/ 77 \ 7T 

min y(k + l)a ± ipo, — J = - 
so the expression becomes 

S(y,<fio) = erf (v^ sin 

p 

+ erf ^ v / 2ysin(ja +v^o)) 

- erf ( \/ 7 h) sin(ja - p 0 )J (28) 

erf ^i/2y sin(fca + <po)) 

+ erf ^-</ 2 y sin (—ka + v?o)) 

+ erf (^\/ 2 y sin((fc + l)a — ip 0 )^j 
+ erf ^ x /2ysin((l - k)a - tp 0 )^ -(29) 


j =i 


_ y4 (-l) fc - 

2-^/ A 


k —0 


Equation (28) matches the one of Dy and Esguerra m 
reminded in Eq. (29 ). 


Then, if the wedge angle is now n/n with n = 2p, 
k = p — 1, so 


min ((fe + l)a + p 0l |) = | 

/ 77 \ 7r 

min [(k + l)a - ip 0 , 7 ; J = {k + l)a - ip 0 = - - ip 0 

and the survival probability can be rewritten 

S(y,<fo) = erf ^^/ 2 ysm^p 0 ^ 
p -1 

+ y](-l) J erf ^v^sin^'a + vJo)) 

3 =1 

-erf [y2ys\n{ja - <^ 0 )) 

+(-l) p erf (v^) - erf (^/ 2 y cos 
/ 2 \ 3/2 

+ (-lf(-) Vy, 


' e * x 


p+00 


du e - pcosh “ sinh | arctan ( ). (30) 


sin (2pyj 0 ) \ 
sinh (pu) J 


This expression can be numerically checked to match the 
known expression of the survival probability for n = 2 

m 

S(y,tpo) = erf( v / 2j/sinv3o)erf( v / 2ycos(po). (31) 
The next values n = 4 and n = 6 have been checked on 


3. Discussion 


We make here several comments on the different terms 
involved in Eqs. (25) and (26) (summarized in Eq. ©). 
Without loss of generality, we assume that the starting 
point is in the inferior part of the wedge (po ^ ck/ 2). 
The first term of Eq. (27) is the survival probability in an 


infinite half space delimited by an absorbing infinite plane 
[T] for a walk whose starting distance from the plane is 
given by ro sin(min(</?o, tt/2)). This latter corresponds 
to the distance of the starting point of the walk to the 
closest absorbing boundary of the wedge (which is the 
one at ip = 0 because of the choice ipo ^ a/2) in the 
original problem. If the wedge is acute, it is rosin^o- If 
the wedge is obtuse, this distance depends on whether 
the projection of the starting point on the axis ip = 0 is 
on the absorbing edge or not (see Fig. [ 3 ]). In the first case 
(corresponding to ipo ^ 7r/2) , this distance is, like in the 
acute case, rosing- In the second case (for ip 0 ^ tt/2), 
the distance is ro, which is the distance from the starting 
point to the apex of the wedge. 

Then, by comparison with the expressions of Dy and 
Esguerra mum, obtained in the particular case of n/n 
wedge angles, the sum involved in Eq. (251 can be seen as 


a sum over generalized images (sinks and sources), that 
only exists for an acute wedge. 

Finally, the integral term and the term fc + 1 of the sum 
of Eq. (271 are the hallmark of a wedge angle different 
from n/(2p +1) with p an integer. 




Figure 3: Distance between the starting point and the wedge. 
If <po ^ 7t/ 2 (left), the projection of the starting point on the 
wedge is on the absorbing edge and the distance is rosim/Jo- 
If po ^ 7r/2, the projection of the starting point on the wedge 
is the apex and the distance is ro- 


D. First-passage time distribution 

Similar expressions for the first-passage time distribu¬ 
tion are easily obtained from Eq. (271. Knowing that 


the first-passage time distribution [IT] (see Eq. (33)). 




(32) 
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it is found that, for any planar wedge, 

fc+1 

+E(- 1 ) J ' x 

3 = 1 

sin (min \ ja + ipo, e 
- sin (min (ja - <p 0 , ) e - ” 


r’o sin^(min(¥3o,7r/2)) 

e 4£)£ 


rg sin^ (minfja + ^o ,7 t/ 2)) 
4 Df 


(min(jQ: —y Q ,7 t/2)) 


+ — / du ( 1 — ^^7 cosiP g ) sinh — e 4 °* 1 


2tt ,/ 0 V 2 Dt 2 


2 X 


arctanf sin (S (^o + f)) V arcta (g fa - f))' 


V sinh (!S) 


V Sinh (fs) 


Using Laplace’s method, we have to distinguish the 
two following cases: (i) when sin (tt (<po + 7r/2) /a) and 
sin ( 7 r (ip 0 — 7t/2) /a) have the same sign, and (ii) when 
they have opposite signs. In the first case, we approxi¬ 
mate the integral by 


± 


e~ 2y 
\/2r ry' 


(37) 


(33) 


In the second case, 

x e ~ 2v ( 1 1 \ 

4ay y sin (g fa + §)) + sin (g fa - §)) J ’ 

(38) 

In both cases, the leading order of the integral term 
is dominated by all the other terms of the sum of error 
functions. Finally, the short time asymptotics of the sur¬ 
vival probability can be defined using a scale of functions 
based on complementary error functions 


with k = (7r/(2a) — 1/2J and y = rfi/fiDt). 


S{t\y,<fio) ~ 1 

y—t 00 


■i>i(y,v 0) 


IV. ASYMPTOTIC DEVELOPMENT OF THE 
SURVIVAL PROBABILITY AT SHORT TIMES 

As an application of the previous results, we now show 
that the asymptotic development of the survival probabil¬ 
ity at short times (y 1) can be conveniently extracted 
from Eqs. (25) and (261. The leading order is 1, as ex¬ 


pected because a walker starting from the bulk cannot 
be on the absorbing boundaries at time t = 0. We are 
interested in the corrections induced by the boundaries, 
and first address the case of an acute wedge. 

Sine being a growing function in [0,7r/2], 

sin(/?o ^ sin (a— ipo) ^ sin(a+^o) ^ sin(2a—<^o) ^ ••• ^ 1- 

(34) 

Moreover, the error function asymptotically grows like 


erf (x) ~ 1 - e~ xZ P ( - 

x —>■ 00 \ X 


(35) 


with P a polynomial. The term erf [y/2ys\n <p$) of 


Eq. (251 thus contains the leading order of the survival 


probability and a first correction to this value, which 
turns out to be the main one. The successive terms of the 


sum involved in Eq. (251 can be shown to be smaller and 


smaller corrections by using the inequalities (341 and the 


asymptotic expansion (351. Last, we evaluate the large y 


asymptotics of the integral term 


X = l — 
7 r 


3/2 


=-y 


/*+ OO 


^~Y 


d u e - y cosh “ sinh — x 

2 


arctan ( sin ^fa + i))\ + arctai Y sin (S G* - f E 


\ sinh(fg) 


V sinh (fg) 


+ Et i’ijjy, Vo) + i’ 2 j+i(y, i Po)\ + R- 2 k+ 2 (,y,vo) 
3 =1 

(39) 

with for j ^ k 

ipi(y,Vo) = - erfc (a/27/sui^o) 

i> 2 j(y, Vo) = (-1) J erfc U/2y sm(ja - <p 0 )) 

ip 2 j+i(y, Vo) = (-1) J+1 erfc ^y/2ysm(ja + Vo)j 

(40) 

and the remainder 
R2k+2(y,<Po) = 

(—l) fc+1 erfc (-faysin (min ((fc + 1 )a — ipo, 


— erfc 


,\ 3/2 


2 y sin (^min ( (k - 

+OO 


1)< 


■ Vo, 


*-y 


+ \n, 


~r 


du e _ycosh “ sinh — 

2 


arctan 


'sin(gfa + §)) 


+ arctan 


sinh(fg) J 

sin(gfa^f)) 
sinh(fg) 


(41) 


and 


( 36 ) 


such that tpi{y,<Po) = ofa-i(y,Vo)), 

R2k+2{y, vo) = o fak+i(y, vo))- 

The obtuse case is simpler. If (pg ^ 7r/2, the error func¬ 
tion term gives the main correction and the integral term 
is subdominant, whereas if (po ^ 7t/2, these two terms 
have the same exponential decay rate. In Fig. [5] we only 
illustrate the first case. 
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Figure 4: Survival probability in a wedge of acute top an¬ 
gle a = 0.4 rad, with y>o =0.1 rad (in red), and the short 
time development cut at different orders. The exact expres¬ 
sion contains here 8 terms, including the remainder Rs(y, <po)- 
We can see that as soon as we keep 2 terms, the short time 
development has a very satisfying range of validity, that can 
be extended by keeping more terms in the development. 



y 


— 1 term 

-exact expression 


V. CONCLUSION 

In this paper, we established simple expressions of the 
survival probability and the first-passage time distribu¬ 
tion in a planar wedge with infinite absorbing bound¬ 
aries. The result holds for any top angle of the wedge, 
and in particular covers the case of obtuse wedges. It 
thus generalizes the expressions obtained by Dy and Es- 
guerra mm, which were limited to wedge angles of the 
form 7r /n with n a positive integer. The final expression 
only involves a finite sum of error functions, that can be 
seen as a sum over generalized images, and an integral of 
elementary functions. 

The expression given here naturally displays a develop¬ 
ment of the survival probability at short times, whereas 
the standard form of the survival probability, that is writ¬ 
ten as an infinite sum of special functions, does not allow 
to get this expansion. Moreover, this short time devel¬ 
opment has a large range of validity. 

The case of biased diffusion in an arbitrary wedge (con¬ 
sidered in m for specific angles of the form n/n) would 
be a natural extension of the formalism developed in this 
work. 

OB acknowledges support from European Research 
Council starting Grant No. FPTOpt-277998. SNM 
acknowledges support by ANR grant 2011-BS04-013-01 
WALKMAT. 


Appendix A: Calculation of the integrals involved in 
the term Ai 


Figure 5: Survival probability in an obtuse wedge in the 
case po ^ tt/ 2 (in red) and the term 1 — erfc (\/2y sin po) 
(in green). The top angle is a = 3.5 rad and ipo = 1.5 rad. 
The error function term provides an accurate approximation 
of the survival probability with a large range of validity. 


The previous analysis shows that at short times (large 
y), for acute wedge angles and obtuse ones where 
tp 0 ^ tt/ 2, the survival probability is mainly influenced 
by the edge closest to the starting point, producing the 
first correction to the limit value 1 

S{y,<Po) ^ 1 “ erfc (\/22/ sin tpo'j . (42) 


We give here details of the calculation of the term Ai , 
which is the sum of the four integrals Cy, + C^~, 

and C£. The first one is 


Ci =2 d ip e 


y cos ?/> 


COS 




dx' 


x ^ cos((2m + l)x') cos ( (2m + l)^- \ (AI) 


m —0 


and using Eq. (151 of the main text, 


Ci = 2 dip e 


y cos 'ip 


ip 

cos- 


(A2) 


Moreover, for acute wedges, our approach also gives a set 
of smaller and smaller corrections, the least correction be¬ 
ing given by the remainder i? 2 fc+ 2 (y, ‘Po)- We check on 
Figs. [4] and [5] that the short time development is accurate 
and has a significant range of validity that increases, in 
the case of acute angles, with the number of correcting 
terms taken into account. In practice, it means that un¬ 
less we need to describe the very long times (small y ), the 
integral term, which is the most complicated to compute, 
can be dropped. 


The integral over x' is 7 t/ 4 if nip/(2a) £ [0, inp^/a}, be. 
if ip £ [0,2iyJo], and 0 otherwise. Thus, as we have by 
definition 2 tp 0 ^ a, the part ip £ [2</? 0 ,ck] gives 0 and 

Cl = f /^ 0#e " C ° S ^ COS 2- (A3) 

We change the variables u = y/2y sm(ip /2) and get 

/7T\3/2 pV _ 

Cl = (-J — erf( v /2ysin^o). (A4) 













by 


The second integral C 2 is the sum of C 2 and C 2 , given As tpo ^ 0, we notice that 


C? =2 Y, 


,y cos(2jct±ip) 


cos [ja ± — 


3=i J0 


„ VQ max( 2 y) 0 ) max (0, (2fc+2)a—7r)) = max( 2 y> 0 ) (2fc+2)a—7r). 

^ da; ' . . ( A9 ) 

Moreover, as — max(a, b ) = min(— a, —b ), we obtain 


+ OO , . V 

x ^ cos((2 m + l)x') cos ( (2 m + l)j7r ± (2 to + 1)^— J Q~ = /ZI\ 

m =0 ' ' ' 2/ 


m=0 
k 

= 2^)(-l ) J I dfi> e 
3=1 Jo 

*< P0 +00 


7 r \ 3 1 4 5 / 2 e y 

7 ^ 


y cos(2ja±i/j) 


COS ja ± 




erf ^y/2ysin ^min ^(fc + l)c 


x / da/ cost(2ro + l)x') cos ( (2m + 1)4— 

Jo t h V '2a 


TYlp 


(A5) 


— erf ^-y/2?/sin (min y(k + l)a — </?o, — j jj (A10) 


We proceed similarly for the last integral: 


Following the same lines, we obtain 


f»max(0,7r— (2fc+2)o;) 


C 3 + = (-l) fc+1 77 


dip e 


y cos((2fc+2)ai+'i/>) 


C 2 = C 2 + c+ = ^(-l) y 


3=1 


7T \ 3/2 e y 

2 > Vy 


ih 

x cos ( (k + l)a + — 




erf(\/2y sin(ja + ip 0 )) - erf(^2y sin(ja - <p 0 )) 

(A6) 


Finally, we compute carefully the last two integrals, 
because of their integration limits 


dx' 5 (x' - (All) 

which is 0 if max(0,7r — (2k + 2)a) = 0. We can then 
rewrite 


/»7T— {2k-\-2)ot 

C+ = (-l) k+1 - / d^e 

z J min(0,7r— (2k-\-2)ot) 


y cos((2fc+2)ai+'i/>) 


cr = (- i) fe+1 A 


di/; e 


2 / cos((2/c+2)a— 1 />) 


x cos ( (jfc + 1 )a + — 


min(0,7r—(2fc+2)o;) 

\ - 77 vo 

v> 


/o 


dx' 6 ( x' — -7- )(A12) 


' max(0,(2fc+2)o;—7r) 
, n( PQ 


x cos I (k + 1 )a — 




dx 6 lx — 


Trip 
2a 


(A7) 


As previously, the integral over x equals to 1 if 
ip £ [0, 2ipo], and 0 otherwise. If the inferior limit is larger 
than 2tpo, the integral is 0. We can then rewrite it as 


7 r r 

C^ = (-l) k+1 ^ / 

z J rr 


max( 2(/?0 5 max ( 0 )( 2 fc+ 2 )a—7r)) 


dV’ 


and end the calculation as before to get 

x erf (i/ 2 ysin (min ((fc + l)a + ipo, 7 ^^ 

— erf (i/ 2 ysin (min ((fc + l)a, ■ (A13) 


max(0,(2fc+2)a—7r) 


/ ? /,\ ouimmng rxjs. 

x e ycos(( 2 k+ 2 ) a -i;) cos I (jfe + i) a _ OL j . (A8) Eq. ([18]) of the main text. 


Summing Eqs. (All, (A6l, (AlOl and (A13l leads to 
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